Multiplets Matter: The Electronic Structure of Rare-Earth Semiconductors and 

Semimetals 
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We demonstrate that a theoretical framework fully incorporating intra-atomic correlations and 
multiplet structure of the localized 4/ states is required in order to capture the essential physics 
of rare-earth semiconductors and semimetals. We focus in particular on the rare-earth semimetal 
erbium arsenide (ErAs), for which effective one-electron approaches fail to provide a consistent 
picture of both high and low-energy electronic states. We treat the many-body states of the Er 
4/ shell within an atomic approximation in the framework of dynamical mean-field theory. Our 
results for the magnetic-field dependence of the 4/ local moment, the infiuence of multiplets on the 
photoemission spectrum, and the exchange splitting of the Fermi surface pockets as measured from 
Shubnikov-de Haas oscillations, are found to be in good agreement with experimental results. 

PACS numbers: 



The established framework for understanding the elec- 
tronic structure of many materials is based on indepen- 
dent electrons subject to a self-consistent potential. This 
description often fails, however, when electronic states 
associated with localized orbitals are involved. Those 
states, while being too close in energy to the valence to be 
considered as 'core' states, do keep a strong atomic-like 
character. For such materials, an entirely different point 
of view is required, in which intra-atomic correlations and 
the full multiplet structure are taken into account from 
the start. Here, we report on a particularly dramatic 
example: erbium arsenide, a semi-metal involving a par- 
tially filled 4/ shell. Interest in this compound was ini- 
tially stimulated by growth, through molecular beam epi- 
taxy, of ErAs onto pseudomorphically compatible III-V 
semiconductor substrates of zinc-blende structure. The 
result is a high-quality, epitaxial metallic contact with 
a continuous anion sublattice through the interfacellf. 
Other growth modes[l|,l2, SlJI have been found to yield 
metallic ErAs nanoparticles epitaxially embedded in the 
semiconductor matrix, with potential applications in- 
cluding efficient thermoelectrics[5| and solid-state THz 
emitters [6j. 

The low-energy electronic states of erbium arsenide 
(ErAs) consist 'li| of a conduction band and a valence 
band, with predominant Er Sd and As Ap character re- 
spectively. The small energy-overlap of these bands is 
responsible for the semimetallic nature of this material. 
In contrast, the Er 4/ electrons are well localized and 
form local moments. Treating simultaneously these two 
kinds of electronic states is a major challenge for con- 
ventional electronic structure theories. Indeed, when 
the 4/ states are treated as core electrons, calculations 
based on the local-density approximation (LDA) do yield 
a semi-metal, but do not otherwise provide an accu- 
rate description of the electronic properties. The vol- 
ume of the calculated FS pockets is approximately three 



times too large, leading to incorrect carrier concentra- 
tions and SdH frequenci es M , Isj compared with accurate 
measurements [qI, 1io|, luJ, ll2] . 

Here we show that a proper treatment of the strong 
correlations and multiplet structure of the 4/ shell pro- 
vides a solution to these difficulties. Our finding is gen- 
eral and has fundamental consequences for the under- 
standing of the electronic structure of a large class of 
rare-earth compounds. ErAs is a particularly remarkable 
example, however, in which the electronic states near the 
Fermi level act as a sensitive probe of the atomic physics 
associated with more strongly bound localized states. 

Our approach is based on the dynamical mean-field 
theory 13] (DMFT), combined[il ^ with LDA. The 
4/ shell is treated as that of an effective atom self- 
consistently coupled to an environment describing the 
rest of the solid. From the hamiltonian of this effec- 
tive atom (which takes into account crystal-field effects, 
intra-atomic Coulomb interactions and the spin-orbit 
coupling), a many-body self-e nerg y is computed within 
the Hubbard-I approximation [16| and inserted into the 
Green's function of the full solid. Self-consistency over 
the total charge density and the effective atom parame- 
ters is implemented. The full description of our approach 
can be found in Ref . |17l | . In the present work it was gen- 



eralized to include the spin-orbit interaction as well as the 
full 4-index local Coulomb interaction matrix. The pa- 
rameter U = 7.94 eV of the local Coulomb interaction on 
the Er 4/ shell has been determined by constrained LDA 
calculations, while the Slater integrals F^ — 12.1 eV, 
_F* = 8.4 eV, and F^ = 6.7 eV, which are known to be 
weakly dependent on the crystalline environment, have 
been taken from the optical measurements on Er ions 
embedded in a LaFa host [l8|. To perform LDA-I-U cal- 
culations we have used the same framework, including 
the spin-orbit interaction, with the Hubbard-I self-energy 
being substituted by a local self-energy computed within 
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FIG. 1: (Color online) Calculated DMFT total magnetic mo- 
ment on Er 4/ shell vs. applied magnetic field at T = 4.2 K 
with (solid blue curve) and without (red dashed curve) crystal 
field effects taken into account. The left inset shows the zoom 
in the field's range up to 20 T. The crystal field splitting of 
the */i5/2 multiplet obtained by LDA+DMFT in paramag- 
netic ErAs is shown in the right inset. 



the rotationally-invariant Hartree-Fock expression. All 
our calculations have been carried out in the rock-salt 
structure at the ErAs experimental lattice parameter of 
5.74 A. 

The ground-state of an isolated Er'^^ ion (Af^^bdP) is 
a 16-fold degenerate multiplet ^/i5/2 {S = §, i = 6, J = 
-^). In a cubic crystal-field and at zero magnetic field, 
this multiplet is expected to split into twofold degenerate 
Fg and F7 multiplets, as well as three 4-fold degenerate Fg 
multiplets 19]. At self-consistency, the zero- field eigen- 
states of our effective atom hamiltonian (inset of Fig. [T|) 
follow this expectation, the obtained ground-state being 
the F7 multiplet. As a magnetic field H is turned on, the 
ground-state remains approximately F7 for a wide range 
of field, and can be decomposed on eigenstates of Jz ac- 



cording to: cii^ L\Jz 



(H) 



Uz 



13\ A- n^^'\ 7" - 5\ 
"13/21''^ ~ 2 I ^ "5/2 \'Jz — 2/ ^" "'-3/2I 

— |) -I- a'_ii/2\Jz — ^'y) ' where the coefficients aj^ de- 
pend on magnetic field. The magnetic moment vs. field 
curve (Fig. [1]) displays an initial sharp rise due to the 
Zeeman splitting of the F7 state, followed by a rather 
slow increase in the range from 5 to 120T due to the pro- 
gressive polarization of the r7 state. From magnetore- 
sistance experiments [12], it was proposed that the mag- 
netic moment is saturated at a value [i20] ~ 5.3/is in fields 
above lOT. Indeed, we find values ranging from 4.8/iB 
at i/ = 5T to 5.6/iB a.t H = lOT, in good agreement 
with the experimental estimate[20|. However, our analy- 
sis shows that this is rather a quasi-saturation, with the 
moment being "frozen" by the crystal field. Eventually, 
at a very high field H ~ 120T we predict a sharp increase 
of the moment to 8.9/iB, corresponding to a symmetry- 
changing transition of the ground-state from r7 to Fg. 

It is interesting to contrast these findings to those of 



LDA-I-U calculations. We performed these for ferromag- 
netic ordering since, for fields above ~ IT relevant to SdH 
experiments, this is the observed configuration [20l| . In 



the LDA-I-U approach [21|, Coulomb interaction effects 



within the 4/-shell are described by a self-consistent, 
orbital- and spin- dependent, one-electron potential. As 
a result, the LDA-I-U solution of lowest energy is found 
to correspond to the filling of independent electron lev- 
els according to Hund's rules [22], leading to the high- 
est possible spin and orbital moments consistent with 
a 4/^^ shell. Indeed, our LDA-I-U results yield a spin, 
orbital and total moment of 2.9fiB, 6.0/xb, and 8.9fiB 
respectively, in contradiction with the observed exper- 
imental value quoted above. Furthermore, as pointed 
out by Larson et al. for nitrides [23], the LDA-f-U so- 
lution breaks the cubic symmetry of the lattice, because 
it is constructed by occupying one-electron eigenstates 
of the angular momentum operator L^. LDA-I-U solu- 
tions preserving cubic symmetry are obtained by occupy- 
ing symmetry-adapted one-electron states (which are not 
eigenstates of L^), but they have higher energy. Hence, 
in a one-electron effective description, the system has 
to choose between minimizing the intra-atomic Coulomb 
interaction energy and preserving the lattice symmetry. 
In reality however, this dilemma does not apply since 
atomic multiplets satisfy both requirements. A proper 
many-body treatment of atomic correlations avoids this 
conundrum, as demonstrated above. We now describe 
the electronic structure of ErAs, obtained with LDA-I-U 
and LDA+DMFT (Fig. [2|). Distinctive features corre- 
sponding to the 4/ states are observed at high energy, 
corresponding to the upper- and lower- Hubbard bands 
(UHB/LHB). In LDA+DMFT, those "bands" truly are 
many-body atomic-like excitations associated with the 
removal (LHB) or addition (UHB) of an electron in 
the 4/ shell. The LHB spans an energy range fronr 
~ — 11 eV to —5.5 eV, in agreement with photoemis- 
sion experiments [H, [11]. Within LDA+DMFT the UHB 
consists not only of the main peak at 2 — 3 eV above the 
Fermi energy, but also of additional multiplet peaks at 
~ 3 — 4eV. As also discussed for (5-PUj25,], this multi- 
plet structure is responsible for the apparent width of 
the UHB, which we find to be ^ 1.5 eV, in agreement 
with inverse photoemmission measurements [23, 24]. 

In contrast, the UHB found with LDA+U spans an 
energy range of only ^ 0.75 eV, because the additional 
multiplet peaks are missed in this approach. Further- 
more, a distinctive feature found within LDA+U is that 
the UHB is almost fully (minority) spin-polarized. This 
is clearly due to the maximal (Hund's rule) spin mo- 
ment, corresponding to the complete filling of the 4/ 
majority-spin states so that no electron addition is possi- 
ble in the majority channel. In contrast, the UHB found 
within LDA+DMFT is found to have only partial spin 
polarization. Hence UHB transitions are possible in both 
spin channels, a distinctive prediction for possible spin- 




FIG. 2: Band structure and density of states (DOS) of ErAs: from LDA+U (left), and LDA+DMFT at an applied field of 5T 
(right). In the corresponding DOS plots the total, Er 5d, Er 4/ and As 4p states are displayed by the black, red, green, and 
blue curves respectively. 



TABLE I: Shubnikov-de Haas(SdH) frequencies (in tesla)and 
carrier concentrations of electrons/holes n (in lO'^^cm"'^). 
The frequencies fsh, fih, and fhh correspond to cross sec- 
tions of the small, light, and heavy hole pockets near the F 
high-symmetry point, respectively. The frequencies e ± and 
e II are due to the transverse and longitudinal cross sections 
of the ellipsoidal electronic pocket. The experimental carrier 
concentrations and SdH frequencies are for an Ero.57Sco.43As 
alloy [12]. 
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polarized photoemission experiments. 

We now show that these differences in the description 
of the high-energy 4/ states have key consequences for 
the electronic properties of states near the Fermi level. 
Near the F high-symmetry point, the As Ap bands form 
"heavy-" {hh) and "Hght-" (Ih) hole pockets, as well as a 
small hole pocket (sh) due to the spin-orbit spHtting of 
the Ap states, while the ellipsoidal electronic pocket at the 
X-point is associated with the Er 5d- states (Fig. [3]). As 
mentioned above, a conventional LDA approach treating 
the 4/ shell as core states severely overestimates the p-d 
overlap, and hence the carrier concentration and the size 
of the FS pockets. Both the LDA-I-DMFT and LDA+U 



TABLE II; Exchange splittings of Shubnikov-de Haas frequen- 
cies A = ./ ^ — f^ (in tesla) 

Orbits Experiment" LSDA-hU GW'' LDA+DMFT 
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"Experimental data for an Ero.57Sco.43As alloy |12|| . 
''Chantis et ai.fH 



treatment of the 4/ states lead to a drastic reduction 
in the magnitude of the carrier density and SdH fre- 
quencies (Table |l|, and hence to a much better agree- 
ment with experiments. These results demonstrate that 
taking into account interaction effects in the 4/ shell is 
key to a proper description of the valence band structure 
of ErAs. Previous works employed an empirical shift of 
the 5d band in order to correct for the inadequacies of 
LDAjSJ, |26|. Despite the overall reduction in the FS 
pocket volume, the LDA+U (self-consistent one-electron 
potential) and LDA+DMFT (self-consistently embedded 
atom) results differ significantly when it comes to the ex- 
change splitting of these pockets. In Table Ull we display 
the measured splittings [Sl, I 111 . Il2| in SdH experiments 
along with the calculated values from LDA+DMFT, 
LDA+U and a recent spin-polarized GW work[27|. It is 
apparent that LDA+U (and spin-polarized GW) strongly 
overestimates the splitting of all orbits, especially for the 
hole pockets, while a much better agreement with exper- 
iments is obtained from LDA+DMFT. The reasons for 
this success are: i) a proper description of the crystal- 
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FIG. 3: (Color online) The low-energy electronic structure of ErAs within the LDA-I-U(top) and LDA+DMFT (bottom) 
approaches. The spin polarization is along the z axis, thus the cubic symmetry is lifted, and the band structure along x(y) 
and z axis are not equivalent. Therefore we designate the corresponding X high-symmetry points as Xx, Xy, and Xz. The 
Brillouin zone is show in insert at the right-hand side of the plot. The overlapping As 4p and Er 5d bands near Ef are shown 
on the left-hand side. They are followed to the right by the cross sections of the hole pockets at the T point in the xy and xz 
planes, and the longitudinal e || (in the TXyW plane) and transverse e _L (in the XyWW plane) cross sections of the electronic 
pocket. The exchange splittings of the pockets are clearly seen. 



field induced freezing of the 4/ spin moment (in contrast 
to the full polarization obtained from LDA-I-U) and ii) 
a proper description of the orbital angular momentum 
content of the upper Hubbard band. Electronic pockets 
are mainly sensitive to the first effect, since the direct 
hybridization between the Er d- and /- states is very 
small near the X-pointjS]. As a result, the splitting of 
the electronic pocket is mainly due to the exchange field 
induced by the polarization of the /-states in the local 
charge density, and hence in the effective potential seen 
by the d-states [28]. Because the spin polarization of the 
/-states is strongly overestimated in LDA-t-U, so is the 
exchange splitting of the electron pockets. In contrast, 
the hybridization between the As p-states and Er /-states 
near the F-point is rather strong. In LDA+U, the UHB 
has a strong spin-polarization in the minority channel: 
hybridization thus repels the corresponding p-states away 
from the Fermi level, leading to too small hole pockets in 
the minority channel (see Table |J). In LDA+DMFT, this 
spin polarization is smaller and the orbital angular mo- 
mentum composition of the UHB is different, with larger 
contributions from components which do not hybridize 
with p-states for symmetry reasons. Hence the smaller 
exchange splitting of the hole-pockets, in better agree- 
ment with experiments. It is remarkable that, because 
of the sensitivity of hybridization effects to both the spin 
and orbital angular momentum of the 4/ states, experi- 
mental measurements of the exchange splitting of the FS 



pockets indirectly probe the electronic state of the Er 4/ 
shell. 

In conclusion, we have shown that a proper descrip- 
tion of the electronic structure of ErAs requires a concep- 
tual framework which gives a central role to local atomic 
physics and multiplet effects, in contrast to conventional 
frameworks based on one-electron effective descriptions. 
This is of general importance for a wide class of mate- 
rials, including all rare-earth based semiconductors and 
semimetals. 
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